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Abstract. In BTA spectra of the supernova SN2006aj, identified with the X-ray flash 
(XRF) and gamma-ray burst XRF/GRB 060218/SN 2006aj, we detected details interpreted 
as hydrogen hues, which is a sign of stellar- wind envelope around a massive progenitor star of 
the gamma-ray burst. Results of modeling two early spectra obtained with the BTA in 2.55 
and 3.55 days after the explosion of Type Ic supernova SN2006aj (z=0.0331) are presented. 
The spectra were modeled in the Sobolev approximation with the SYNOW code (Branch et 
al. 2001; Elmhamdi et al. 2006). In the spectra of the optical afterglow of the X-ray flash 
XRF/GRB 060218 we detected spectral features interpreted as (1) the Ha PCyg profile for 
the velocity 33000 km s~^ — a wide and almost unnoticeable deformation of continuum in 
the range of ~ 5600 — 6600A for the first epoch (2.55 days) and (2) a part ("remnant") of 
the Ha PCyg profile in absorption blueshifted by 24000 km s^^ — a wide spectral feature 
with a minimum at ~ 6100A (the rest wavelength) for the second epoch (3.55 days). Taking 
into consideration early BTA observations and spectra obtained with other telescopes (ESO 
Lick, ESO VLT, NOT) before 2006 Feb. 23 UT, it can be said that we observe evolution of 
optical spectra of Type Ic core-collapse supernova SN 2006aj during transition from the short 
phase related to the shock breakout to outer layers of the stellar- wind envelope to spectra of 
the phase of increasing brightness corresponding to radioactive heating. Signs of hydrogen 
in spectra of the gamma-ray burst afterglow were detected for the first time. 



I. INTRODUCTION 

On Feb. 18.149, 2006 UT the space observa- 
tory Swift detected a peculiar gamma-ray burst 
(GRB) with a powerful component of supernova 
(SN) emission in spectra and in the light curve 
of the GRB afterglow. Therefore, this burst is 



simultaneously classified both as GRB 060218, 
and as SN2006aj: GRB 060218/SN 2006aj. But 
since in this case the X-ray emission was prevail- 
ing in the GRB spectrum, the GRB is also clas- 
sified as XRF (X-Ray Flash), and the event is 
denoted either as XRF/GRB 060218/SN 2006aj, 
or as XRF 060218/SN 2006aj. Below we will 
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mostly use the latter notation to emphasize the 
circumstance which is important for this work 
that in the case of SN2006aj the SN event it- 
self started with observation of a powerful X-ray 
flash (XRF). 

For this XRF/GRB, direct observational 
signs of an early phase of expansion of shock 
arising due to explosion of a compact core of 
a pre-SN and the shock breakout to the outer 
boundary of the stellar-wind envelope surround- 
ing the SN progenitor star were obtained for the 
first time. During the first 2 hours, the shock- 
wave was observed as the mostly X-ray flash 
XRF 060218, and afterwards, when the envelope 
became optically thin, as a powerful ultraviolet 
(UV) burst with the maximum brightness in 11 
hours after GRB detection (Campana et. al., 
2006; Blustin 2007). During first 2800 seconds, 
beside non-thermal emission typical of GRB af- 
terglows, XRF 060218 also demonstrates in its 
X-ray spectrum a powerful thermal component 
whose temperature falls, and with time the emis- 
sion maximum shifts to the UV and optical spec- 
tral range (see Fig.l in Campana et. al., 2006 
and Fig. 2 in Blustin 2007). 

XRF/GRB 060218 is one of the nearest GRBs 
with its redshift 2=0.0331. In this respect it can 
be compared to GRB 030329/SN 2003dh having 
the redshift 2;=0.1683 which was also identified 
with the Type Ic SN. Both events have aroused 
considerable interest because such coincidences 
(GRB and SN) happen rarely, only once every 
two-three years (Chapman et al. 2007), when 



GRBs revealed SNe starting from the moment 
of star core explosion. It can be said that the 
study of GRBs is a new phase of investigation 
of the same core-collapse SNe, but from the very 
beginning of this remarkable phenomenon. That 
is why in these cases early spectral observations 
turn out to be very important for understanding 
the mechanism of both core-collapse SN explo- 
sion itself and GRB source. When such early ob- 
servations are successful, the 6-meter telescope, 
with its large accumulating area and eastern 
location relatively other large European obser- 
vatories, can play an important role in imple- 
mentation of the international program of spec- 
troscopic observations (monitoring) of quickly- 
fading optical GRB afterglows. 

Fig. 1 represents the spectra obtained 
with the 6-meter telescope (Fatkhullin et al. 
2006). Like the case of GRB030329/SN2003dh 
(Sokolov et al. 2003), the spectra of 
XRF060218/SN2006aj are among the earhest 
spectra obtained with a high signal/noise ratio 
(as is shown in Table 1), which permits apply- 
ing the interpretation methods usually used for 
spectra of Type Ic SNe (Branch et al. 2001; 
Elmhamdi et al. 2006). 

As was shown by Campana ct al. (2006), 
both the XRF itself and the UV burst in about 
10-11 hours after XRF 060218, and then the UV 
excess in early spectra of the afterglow can be ex- 
plained by interaction between the SN shock and 
a stellar-wind envelope around the massive pro- 
genitor star. This is a so-called "shock breakout" 
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effect (the shock breaks out through an envelope 
around a collapsing and exploding star core). 

Such an effect for core-collapse SNelb/c and 
SNell has been known for a long time already 
(Colgate, 1968; Imshennik and Nadezhin, 1989; 
Calzavara and Matzner, 2004). It can be ob- 
served as a relatively short phase of the SN ex- 
plosion which, after beginning with an X-ray 
flash (XRF), ends with a bright UV burst, which 
announces arrival of the shock to surface of the 
exploding star. It can be also said about the 
shock breakout to upper layers or to an optically- 
thin "surface" of the extended stellar-wind en- 
velope surrounding the collapsing and explod- 
ing cores of Type lb or Ic SNe. Below, the 
term "shock breakout" will be used. Like two 
famous SN 1993J and SN 1987A, in the case of 
XRF060218/SN2006aj this effect was observed 
as a sharp and short early peak in the optical 
light curve — see Fig. 2 in Campana et al. 2006 
and Fig. 2 in this paper. We mentioned the same 
effect when explaining the very first spectra of 
the GRB030329/SN2003dh afterglow obtained 
in ~10 hours after GRB 030329 (Sokolov 2003). 
Due to detecting XRF060218/SN2006aj in the 
beginning of the SN burst, i.e. before the shock 
breakout, we managed to observe motion of the 
shock inside the stcllar-wind envelope (the shock 
breakout effect) during the first two hours as an 
X-ray fiash (XRF) with thermal spectrum. 

From Table 1 it is seen that BTA spectra 
refer to minimum in the light curve, i.e. we 
observe the transition from phase 1 of thermal 



emission related to the shock breakout to "sur- 
face" of the stellar-wind envelope (the first peak 
in Fig.2) to phase 2 of the subsequent increasing 
of SN brightness with maximum in ~10 days in 
Fig.2, which corresponds to the (non-thermal) 
radioactive heating of the expanding SN enve- 
lope caused by the decay ^^Ni^^^Co— >^^Fe. 

Thus, during our observations, physical con- 
ditions have been changing quickly, and this pa- 
per is dedicated to interpretation of spectra of 
this remarkable transitional phase. But before 
comparing observational and theoretical spectra 
in Section 2 of this paper, we adduce estimates 
of typical energy, temperature, size and velocity 
which directly follow from results of Swift ob- 
servations of XRF 060218 in the X-ray and UV 
ranges obtained before the beginning of the BTA 
spectral observations. 

1) Energy of the X-ray and UV flashes: 
Light curves of the XRT (0.3-10 kcV) and 
UVOT ranges are given in the paper by Cam- 
pana et al. (2006). Ibidem, there is an estimate 
of energy radiated in gamma- and X-rays during 
the first two hours when the shock has been still 
moving (breaking out) through the wind. En- 
ergy release in the X-ray (XRT) flash is equal to 
~ 6 X 10"^^ erg. Then, in about 8-11 hours, the 
UV (UVOT) light curve shows a powerful peak 
(the first maximum in our Fig.2) caused by the 
shock breakout to the surface or, more exactly, to 
the outer boundary of the stellar-wind envelope 
when it becomes sufficiently transparent. It is 
this UV burst that is the Colgate shock break- 
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out effect (Colgate 1968). Energy released in this 
UV flash in the Swift/UVOT range during 28 
hours and estimated from data of the paper by 
Campana et al. (2006) gives the number of the 
same order ~ 3 x 10*^^ erg, which directly says 
about identical nature of the X-ray (XRT) and 
UV (UVOT) flashes/bursts. 

2) Evolution of temperature: 

At first, when the shock breaks out through 
the wind envelope and energy is released 
mostly in X-rays, the thermal spectrum corre- 
sponds well to the temperature kT ~ 0.17 keV 
(~2xl0^ K). The temperature falls to 0.03 - 0.05 
keV (350000 - 580000 K) or lower, by the time 
of the UV (UVOT) flash in 11 hours after the 
GRB (the first maximum in our Fig.2) taking 
into consideration the kT measurement uncer- 
tainties with Swift/UVOT in this period (only a 
part of the black-body peak gets into the UVOT 
range, see Blustin, 2007). Then the tempera- 
ture of the thermal component decreased down 
to -43000 K or to kT 3.7 -M.9/-0.9 eV (Cam- 
pana et al. 2006, Fig. 3) in the subsequent Ri22.5 
hours. Thus, it is thought that the temperature 
can be even lower, less than 10000 K, by the mo- 
ment of our first spectral observation (2.55 days 
after GRB, see Table 1), in accordance with esti- 
mates of the temperature decrease rate adduced 
in the paper by Campana et al.. 

3) Size of the stellar-wind envelope: 

Prom the data given in the paper by Cam- 
pana et al. (2006), one can also estimate ra- 
dius of the wind envelope surrounding the Wolf- 



Rayet (WR) star before the explosion. (For ex- 
ample, the arguments in favor of the WR star 
surrounded by an extended stellar- wind envelope 
as a progenitor star of XRF060218/SN 2006aj 
are given by Blustin, 2007). It is natural to 
connect the size of this envelope to the bright 
UV (UVOT) flash observed in ~11 hour after 
GRB. At this moment the shock which was ob- 
served before that only in X-rays becomes vis- 
ible in UV and optical at last, since layers of 
the wind envelope over the shock become opti- 
cally thin and the shock breaks out to the "sur- 
face" (more exactly, to upper layers) of this ex- 
tended envelope related to the GRB progenitor 
star or pre-SN. At this moment, according to 
Swift/UVOT data about evolution of tempera- 
ture and radius of the thermal component of the 
GRB/XRF 060218 afterglow, i.e. at a temper- 
ature of kT ~ (0.03— 0.05) keV (</~580000K) 
at the moment (~ 10*^ s) of the luminosity max- 
imum, the size of the wind envelope around the 
pre-SN must be equal to >/ ~3OOi?0 at the bolo- 
metric luminosity (4.6 — 35.5)xl0^^ erg/s. 

4) Velocity of the shock: 

While the shock moves inside this stellar- 
wind envelope, the radius corresponding to the 
thermal radiation component (related to the 
shock breaking out from the center) is constantly 
increasing from ^ 5.7 Rq to Ri 17i?o, as it follows 
from the same data (Campana et al., 2006) of the 
X-ray observations carried out with Swift/XBT 
before maximum of the UV flash. The radius of 
the thermal component continues increasing to 
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« 47OOi20 in 1.4 days (see Fig. 3 in Campana et 
al.), i.e. by the time of observation carried out 
with Swift/\J\OT after the flash maximum. At 
this time the shock has already broke out to the 
"surface" of the stcllar-wind envelope and lumi- 
nosity of the thermal component has decreased 
(see Fig.2). Dividing the path of the shock (or 
the radius ^3.29 +0.94/-0.93 x 10^^ cm of the 
thermal component at kT ^ 3.7 +1.9/-0.9 eV) 
by the time (1.4 days) one can obtain expan- 
sion velocity of the photosphere related to the 
shock by the moment of the spectral observa- 
tions: (2.7 +/- 0.8) X lO^cms-i (Blustin 2007). 
Such a velocity is typical of core-collapse SNe, 
and it is comparable with line widths observed 
in the spectrum of SN2006aj (Plan et al. 2006). 

Further it will be shown that the above esti- 
mates of the energy, size and velocity obtained 
from the observations of the shock breakout ef- 
fect by the Swift /XRT /\JYOT observatory are 
also confirmed when analyzing our optical spec- 
tra of XRF060218/SN2006aj obtained in 2.55 
and 3.55 days after the beginning of the SN ex- 
plosion, i.e. when contribution of the thermal 
component of radiation of the shock was still 
determining, which is also indicated by strong 
blue excesses in our spectra (Fig. 1). But, as 
is seen from the UBVR light curves (Fig.2), in 5 
days the GRB afterglow has been already notice- 
ably reddening, which is related to the change of 
the mechanism of the SN radiation by this time, 
when the classical photospheric phase of the SN 
envelope expansion begins. This phase is usually 



observed in the core-collapse SNe and it was de- 
scribed in reviews many times (see, for example, 
Imshennik and Nadezhin, 1988). 

Section 2 describes early spectral observa- 
tions of the XRF 06021 8/ SN 2006aj optical after- 
glow carried out with the 6-meter SAO RAS tele- 
scope. Section 3 describes the modelling of spec- 
tra with the SYNOW code and discusses spec- 
tral manifestations of the wind envelope as wide 
absorptions near 5800 A and 6 100 A (rest wave- 
lengths with z=0) corresponding to the PCyg 
absorption in the Ha hydrogen line. Section 4 
contains the discussion of results in the context 
of: 1) an obvious evolution of the Ha traces by 
data for XRF060218/SN 2006aj obtained with 
BTA and other telescopes; 2) manifestations of 
the Colgate shock-breakout effect in the light 
curves, spectra, sizes and luminosities of "usual" 
Type Ib/c SNe; and 3) asymmetry of explosion 
of the core-collapse SNe. Conclusions are formu- 
lated in Section 5. 

II. EARLY SPECTRAL OBSERVATIONS 
OF THE OPTICAL AFTERGLOW OF 
XRF060218/SN2006AJ WITH THE 
6-METER SAO RAS TELESCOPE 

Within the context of the program of follow- 
up observations of GRB optical afterglows, 
the variable object detected with Swift/UVOT 
and associated with the XRF 060218 event was 
spectroscopically observed with the 6-meter 
SAO RAS telescope. We used the Spectral 
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Camera with Optical Reducer for Photo- 
metrical and Interferometrical Observations 
(SCORPIO) mounted in the BTA primary 
focus [http : / / www . sao. ru/hq/lsfvo / devices / 
Scorpio/ scorpio.html]. There were two sets of 
observations: at the epochs 16:31-17:18, Feb. 
20, 2006 UT and 16:31-17:16, Feb. 21, 2006 UT. 
The dispersing element was a combination of a 
transparent grid and the prism VPHG550G with 
the spectral operational range and resolution 
(FWHM - Full Width at the Half Maximum) of 
3500 - 7500AA and lOA, respectively. 

The processing of data was standard. It in- 
cluded subtraction of electronic zero, correction 
for flat field, wavelength calibration with the 
help of comparison spectrum of a Ne-Ar lamp, 
correction for atmosphere extinction and cali- 
bration by an absolute flux with the use of ob- 
servations of a spectrophotometric standard at 
every night. To detect a possible short-time 
variability and to control calibration of absolute 
fluxes from night to night, a bright stellar- like 
object was set at the input slit together with 
the optical transient (see Fig. [3|). In all, 4 expo- 
sures were obtained in the first and second peri- 
ods. We did not detect any short-time variability 
within these observational periods, therefore we 
summed all 4 individual spectra to get a higher 
signal/noise ratio. Figure [U shows average spec- 
tra in every night, Table [H lists average epochs 
of these spectra. 

Then observational spectra (Fig. 1) were cor- 
rected (Figs. 5 and 6) for galactic extinction ac- 



cording to the dust distribution maps of Schlegel, 
Finkbeiner, & Davis (1998). The value of 
E{B — V) in the direction of the object (a2ooo = 
03'^21'^39^683, (52ooo = -hl6°52'01".82) was 
equal to 0.14. When accounting for extinc- 
tion, the dust-screen model was accepted, in 
which the expression for extinction has the form 
F^ntiX) = F,,,(A)100-4 '=W^{^-^), where F,nt{X) 
and Fohs{X) are issued (without the extinction) 
and observed fluxes respectively. The extinction 
curve of the Milky Way {k{X)) was taken from 
the paper Cardelh, Clayton, & Mathis (1989). 

III. THE COMPARISON OF THE 
OBSERVATIONAL SPECTRA OF 
XRF060218/SN2006AJ WITH 
SYNTHETIC SPECTRA 

To interpret spectra we used the "SYNOW" 
code calculating synthetic spectra of SNe 
(Branch et al. 2001). This multiparametric 
code was already applied many times for a di- 
rect analysis of spectra of core-collapse SNe, 
when hydrogen was detected in the phase of 
radioactive heating of the expanding SN enve- 
lope (Branch et al. 2002, Baron et al. 2005, 
Elmhamdi et al. 2006). The computation algo- 
rithm of model spectra is based on the following 
assumptions: spherical symmetry, a homologi- 
cal expansion of layers and a sharp photosphere. 
When interpreting early spectra, this optically- 
thick photosphere emitting a continuum black- 
body spectrum in some cases can be identified 
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with the shock. At least this is vahd in the first 
days (Tsp=1.95d and 2.55d from Table 1) after 
XRF 060218, when spectra (Fig.l) and photome- 
try (Fig. 2) still demonstrate a large blue excess, 
which is related with influence of the thermal 
component of radiation from the first short and 
powerful UV flash. It is assumed that spectral 
lines are formed over this expanding photosphere 
as a result of resonance scattering which is inter- 
preted in the SYNOW code in the Sobolev ap- 
proximation (Sobolev 1958, Sobolev 1960, Gor- 
batsky and Minin 1963, Sobolev 1967). 

This code allows us to identify lines, to esti- 
mate velocity of the photosphere expansion and 
the velocity range (or scatter) for lines of ev- 
ery ion detected in the spectrum. The detailed 
description of the SYNOW code can be found 
in the fore-quoted papers Branch et al. 2001, 
Elmhamdi et al. 2006. To understand our re- 
sults, the fact is essential that in the SYNOW 
code the velocity of layers expansion over the 
photosphere is proportional to distance from ev- 
ery point of an expanding layer to the center 
{v ~ r). 

Judging by the form of observed spectral 
features, there are two cases: the undetaching 
( "the undetached case" ) and detaching ( "the de- 
tached case") of the layer, in which a spectral 
line forms, from surface of the expanding sharp 
photosphere. If matter of moving gas layers lo- 
cated over the photosphere does not detach from 
it (i.e. all layers radiate and screen photosphere 
starting from its level) , then a line is observed in 



absorption or in emission like the case of here- 
considered case (Fig. H]) of a line of the classic 
PCyg profile. We have chosen this model of 
layers undetached from photosphere ("the un- 
detached case") to interpret our first spectrum 
(Tsp=2.55d, see Fig. 1). Besides, apparently, at 
this moment the expanding photosphere still can 
be identified with the shock, because our first 
spectrum (with a larger excess in its blue part) 
was obtained closer to the moment when the 
shock breaks out to the "surface" of the stellar- 
wind envelope surrounding the pre-SN — the 
first maximum in Fig. [2j 

A layer detached from the photosphere man- 
ifests itself in spectrum as a smoothed emis- 
sion and a strongly blue-shifted absorption — 
a "remnant" (part) of the PCyg profile — as is 
shown in Fig. |4l The "detached case" is realized 
when velocity of the layer in which the line forms 
noticeably exceeds velocity of the photosphere. 
If t> ~ r, then the absorption part of the PCyg 
profile is a result of screening photosphere by a 
narrow layer; and the narrower is this layer (radi- 
ating a spectral line), the closer is the emission to 
the photosphere continuum. Such a model cor- 
responds to our second spectrum (Tsp=3.55d), 
at least, for HI lines which can be connected to 
the wind envelope. It may be assumed that the 
shock provoked the motion and "detachment" of 
the uppermost layers of this extended wind enve- 
lope surrounding the GRB/SN progenitor star. 

In 3.55 days after the SN explosion, when the 
contribution of the thermal component deceases 
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(Fig. 2) and the expanding envelope of SN be- 
comes more and more transparent, the shock 
cannot be already identified with the photo- 
sphere totally. The so-called photosphere phase 
of the SN explosion begins. The SYNOW code 
was applied mostly to interpret spectra of this 
phase (Elmhamdi et al. 2006, see also references 
ibidem) . 

In Figs. [5] and [6] the BTA spectra of 
XRF060218/SN2006aj, obtained in 2.55 and 
3.55 days after the event are compared to syn- 
thetic spectra modeled with the SYNOW code. 
The first spectrum taken right after the bright 
UV burst (Fig. [2]), is easily modeled for a rather 
simple set of parameters. We have chosen tem- 
perature of the photosphere (tbi,) used in fit- 
ting the energy distribution in the spectrum in 
Fig. [5] (rest wavelengths) to be equal to 9000 K. 
This corresponds also to the temperature fall 
rate by results of Swift /XRT/UVOT observa- 
tions (Campana et al. 2006), according to 
which by the time of our spectral observation 
the temperature must be less than 10000 K. In 
the spectrum taken in 2.55 days the SYNOW 
model shows the photosphere velocity equal to 
33000 km s^^. This parameter is within errors 
of estimates of the photosphere expansion veloc- 
ity connected with the shock, which is obtained 
from the same Swift /XRT/UVOT data ((2.7 +/- 
0.8)xl0^ cms"^). 

This means that by the moment of our spec- 
tral observations and approximately in a day af- 
ter the last Swift /UVOT observation of the fad- 



ing UV flash (the shock breakout), the shock ve- 
locity remained within the same limits. Further- 
more, a wide and hardly noticeable depression 
in continuum (see Fig. [1] with observed wave- 
lengths) at ~ 5900 — 6300^, and an almost un- 
noticeable excess of flux within the wavelength 
range f« 6300 — 6900A is best fitted by a wide 
Ha PCyg profile (see Fig. [5]) at the same ve- 
locity 33000 kms~^ In the SYNOW code this 
corresponds to "the undetached case" in Fig.dl 

It may be supposed that at this time a part of 
the stellar-wind envelope located over the pho- 
tosphere or an extended layer over it moves to- 
gether with the photosphere (the shock) , and we 
observe acceleration or detachment of upper lay- 
ers of this envelope. This is one more result (in 
~2 days after the bright UV (UVOT) flash) of 
the shock breakout to "the surface" of the mas- 
sive circumstellar envelope which was (almost) 
resting before the SN explosion. According to 
estimates given in Introduction, the size of this 
envelope was >/~300i?o, at least at the mo- 
ment of the shock breakout in ~11 hours after 
the event. To demonstrate that an almost un- 
noticeable oscillation of flux in Fig. [1] within the 
range ~ 5900 — 6900A of the observed spectrum 
may be indeed a very wide Ha PCyg profile for 
velocities about 0.3x10^*^ cms~^, in Fig.[5]we give 
a model spectrum with the classic PCyg profile, 
but for a much less velocity equal to 8000 km 
s^^. (The laboratory Ha wavelength is marked 
by a narrow emission of the host galaxy.) 

We made fitting of this undetached case by 
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different syntlietic spectra, but velocity of the 
photosphere (Vphot), all elements and their ions 
(vmine) was identical and equal to 33000 km 
s~^. Basic parameters of calculation are given 
in Table 2 for the "undetached case" in Fig. [5l 
The only difference is in effective depths of ion 
line formation r (taul), which can be chosen 
relatively small (r < 0.5) for all elements and 
their ions, for the observed spectrum to be de- 
scribed satisfactorily. Within the wavelength 
range ~ 4500 — 7200 A it is hydrogen that has the 
largest value r=0.2. Other elements and their 
ions were taken into account in the synthetic 
spectrum with even smaller r, which may reflect 
small variations in relative abundance of other 
elements in comparison with the most abundant 
one. This may be connected with the fact that 
in this first spectrum {Tsp=2.55 days) the con- 
tribution of the stellar-wind envelope is still de- 
terminant. 

It should be said that we were trying to get 
the best agreement between the observed and 
synthetic spectra in their range where observa- 
tions give the best signal/noise ratio and the 
theoretic spectrum is within the grass or in its 
nearest proximity. In particular, the outlining 
deficit of flux to the left of the emission line 
of the host galaxy [Oil] 37271 (in Fig.[5]) can 
be described by influence of Call with r=0.5. 
This is the largest values of the effective depth 
of formation that we used in this paper dedi- 
cated to interpretation of the earliest spectra of 
the XRF060218/SN2006aj afterglow. 



In our second spectrum for Tsp=3.55 days 
(see Fig.[T] with observed wavelengths) we see 
some absorption at 6300A interpreted here as a 
strongly blue-shifted "remnant" of the Ha PCyg 
profile. The best fitting of the synthetic spec- 
trum to the observed one corresponds to the case 
denoted as "the detached case" in Fig. HI when a 
part of the wind envelope has already detached 
from the photosphere. I.e. contribution of the 
thermal burst (Fig. 2) is quickly decreasing and 
the shock is becoming more and more transpar- 
ent — the photosphere phase of the SN envelope 
expansion begins (Branch et al. 2001). 

We were modelling this second spectrum 
within the velocity range 18000 km < V < 
24000 kms~^ (see Fig.[6]of rest wavelengths with 
z=0). The synthetic spectrum which is the clos- 
est to the observed one is the spectrum with pa- 
rameters of the model from Table 3, where the 
photosphere velocity is equal to 18000 kms^^, 
and its temperature is t^f, = 8200 K. In the syn- 
thetic spectrum we took into consideration the 
contribution of lines with small t(< 0.3) of such 
elements as HI, Hel, Fell, Sill, 01, Cal, Call, 
Till, NI, CI, CII, Mgl, Mgll, Nal. Fig.[6] shows 
location of some of these lines in those spectral 
regions where the contribution of the ion into 
the spectrum is essential with given parameters 
of the model from Table 3. The list of lines used 
in the SYNOW code ("the reference lines") is 
given in Elmhamdi et al. (2006), detailed in- 
formation is given directly in the SYNOW code 
http: / / www.nhn.ou.edu / ~parrent / synow.html. 
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A wide absorption with minimum near 6IOOA 
(Fig. [6]) can be described by influence of HI for 
"the detached case" (Fig.|4]) at r = 0.16. In this 
region (~ 5700 — 6500A in Fig. [6]) the synthetic 
spectrum is characterized by a smoothed emis- 
sion from the red side and a strongly blue-shifted 
absorption with minimum near 6IOOA — "the 
remnant" of the Ha PCyg profile of the line HI. 

Thus, in Tsp=3.55 days, hydrogen has al- 
ready detached from the photosphere and the 
corresponding layer moves at a velocity of 
24000 km s~^ (see Table 3). Here also all ions 
were taken with the small values r < 0.3. But 
we failed to describe totally this observed spec- 
trum with identical values of all other parame- 
ters, as was made in Table 2 ("the undetached 
case" ) . The layers where lines of other elements 
mainly form either move at the velocity identical 
to that of the photosphere (18000 km s^^ for the 
lines of 01, Cal, Call, Till, NI, Nal), or at the 
velocity of hydrogen (24000 km s'^ for the lines 
of HI, Hel, Fell, Sill, CI, CII, Mgl, Mgll). For 
the best description of the observed spectrum, 
typical velocity values (the parameter ve in Ta- 
ble 3) determining the characteristic thickness of 
layers occupied by every element should also be 
chosen different. It may be assumed that this 
spectrum is already affected more by the chem- 
ical composition of the scattering SN envelope 
(or envelopes) which changed result of evo- 
lution and explosion of a far-evolving massive 
core of the progenitor star. 



IV. DISCUSSION OF THE RESULTS 

4.1 Evolution of spectrum of SN 2006aj 
and other core-collapse SNe. 

A wide and low-contrast feature in the spec- 
trum of the XRF060218/SN2006aj afterglow 
with minimum at ~6100A (in rest wavelengths) 
can be traced in all early spectra obtained with 
our and other telescopes starting from Feb 21.70 
UT (see Table 1). The spectrum which is the 
closest in time (Feb 21.93 UT) to our second 
spectrum from Table 1 is the spectrum obtained 
with NOT (see Fig. 2 in Sollerman et al. 2006), 
which also contains minimum at the same wave- 
length ~ 6IOOA. The same feature is also seen in 
the ESO Lick spectrum obtained on Feb 22.159 
UT, but its signal/noise ratio (it is the second 
spectrum in Fig. 1 in Mazzali et. al. 2006) 
is noticeably less than for two above-mentioned 
spectra. From data adduced in Fig.l in the pa- 
per Mazzali et. al. (2006) it is seen that this 
feature obviously evolves starting from the very 
first VLT spectrum and becomes the deepest in 
the spectrum obtained with VLT on Feb 23.026 
UT, i.e. in about 5 days after XRF 060218. 
Later spectra obtained after minimum in the 
light curve (see our Fig. [2]) are more and more 
affected by the increasing wide absorption re- 
lated to Sill near 6OOOA, as was noted by Maz- 
zali et. al. (2006). But the narrow minimum 
at 6100^ is still seen even in the March's VLT 
spectra. This is very similar to what is observed 
in spectra of some Type lb SNe (Parrent et al. 
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2007), where the narrow absorption related to 
Ha is also detected and evolves in the spectra 
obtained in the beginning of the long rise iden- 
tical to that in our Fig. [2] to the typical (of Type 
lb and Ic SNe) maximum of brightness. 

The start of evolution of the spectral 
feature at 6100A which is then traced for 
XRF060218/SN2006aj in the data from other 
telescopes is also confirmed in our above- 
mentioned spectra. But we also include here 
the first BTA spectrum of Feb 20.70 interpret- 
ing the wide and almost unnoticeable depres- 
sion of continuum within the range 5600 — 6600A 
(5800 — 68OOA are observed wavelengths in Fig.[T] 
for Tsp=2.55d) as the Ha PCyg profile for veloc- 
ities of 33000 km s"^ The same small oscillation 
of continuum can be seen also in the very first 
VLT spectrum of Feb. 21.041 UT (see the same 
Fig. 1 in Mazzali et. al. 2006) obtained in 8 
hours after our first spectrum. As to the MDM 
spectrum obtained in 14.5 hours before us (see 
Table 1), here also one can see a weak change 
of continuum at ~ 5700A which is noticeable 
even at the low signal/noise ratio (Mirabel et 
al. 2006). Here also the same very wide Ha 
PCyg profile for velocities of --^30000 km s^^ is 
quite competitive to the identification suggested 
by these authors. Thus, taking into account 
all early observations from Table 1, it may be 
said that we do observe evolution of the optical 
spectrum — a transition from the phase of the 
Colgate shock breakout effect and the powerful 
(thermal) burst related to it to spectra of the 



phase of the SN brightness increasing which cor- 
responds to radioactive (non-thermal) heating at 
the decay ^^Ni^seQ^^sepg^ 

Maybe, this is the main difference and nov- 
elty of our approach to interpretation of early 
spectra of XRF060218/SN 2006aj to that repre- 
sented in the paper by Mazzali et. al. (2006), 
where the VLT/Lick spectra of this SN were 
identified and where the question is only on the- 
oretic spectra for the phase of radioactive heat- 
ing — see the light curve (Fig. 2) in their paper. 
And when analyzing spectra of usual Type Ib-c 
SNe ("stripped-envelope SNe"), the authors of 
the SYNOW code (Branch et al. 2001, Branch 
et al. 2002, Baron et al. 2005, Elmhamdi et 
al. 2006) were calculating their synthetic spec- 
tra with traces of the Ha line, but yet not for so 
early phases with a considerable contribution of 
the thermal radiation from the UV flare and so 
large expansion velocities which were observed 
for XRF 060218/SN 2006aj. 

The evolution of ah 16 ESO Lick and ESO 
VLT spectra of SN2006aj (Pian et. al. 2006) 
was modeled up to March 10, 2006 (20 days 
after XRF 020618) in the fore-quoted paper by 
Mazzali et al. (2006) with the help of a more 
sophisticated code of synthesis of SN spectra 
by the Monte Carlo method based under the 
same assumptions as the SYNOW code, but with 
consideration of model distributions of density 
and temperature in the envelope over the pho- 
tosphere, radiation transfer in terms of transi- 
tions in lines and electronic scatting (details of 
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the method are expounded in the papers Maz- 
zah and Lucy (1993), Lucy (1999), Branch et 
al. (2003). All typical signs of Type Ic SNe in- 
creases in spectra before the wide maximum in 
the SN2006aj light curves (psIO days after GRB 
in Fig. [2]) and after it. The spectra were mod- 
eled for the velocity range 20000 km s"^ < v < 
30000 km s~^. The strongest features mentioned 
by the authors are the lines of Fell, Till, and at 
later phases — Call (< 4500l), Felll Fell (near 
50001), Sill (near 6OOOA), 01 (near 7500l) and 
Call (near 8OOOI). 

We also took into consideration all this list 
of lines when interpreting our early spectra with 
the SYNOW code (see Fig. [6]) for almost identi- 
cal range of velocities, but we also included con- 
tributions of the HI He I lines. From calcula- 
tions adduced in Mazzali et al. (2006) it is seen 
that in the very early spectra of Feb 21, Feb 22, 
Feb 23 the authors did not take into considera- 
tion abvious traces of the Ha line, and in later 
spectra only the increasing influence of Si II near 
6OOOA was considered — see Fig. 1 in Mazzali et 
al. (2006) and the caption to it. At the same 
time, as was mentioned above, the absorption 
with minimum at ~ 6IOOA can be traced in the 
observed spectra up to March, 4. It is the most 
noticeable feature near 6OOOA in the early spec- 
trum of Feb. 23 in the whole range from ~ 5000A 
to ~ 8500^, where the Monte Carlo method does 
not fix any absorption at ~ 6IOOA. 

4.2 The signs of hydrogen in core- 
collapse SNe spectra. 



The signs of hydrogen in spectra of Type 
lb and Ic (Ib-c) SNe are no news. In particu- 
lar, the signs of hydrogen and evolution of the 
blueshofted Ha line were already found with the 
help of the same SYNOW code in the analysis 
of a time series of optical spectra for usual core- 
collapse Type Ic and lb SNe. In this analysis, 
special attention addressed to traces of hydro- 
gen in observations of especially these stripped- 
envelope Type Ib-c SNe. Though, by (formal) 
definition, the Type Ic and lb supernovae does 
not have conspicuous lines of hydrogen in its op- 
tical spectra. The Ib-c SNe are usually modelled 
in terms of the gravitational collapse of massive 
and bare carbon-oxygen cores which stripped en- 
velope before collapse, and, apparently, signs of 
this envelope must be present always in spectra 
of these SNe as hydrogen lines. Though most 
often hydrogen can be identified more or less re- 
liably only in sufficiently early spectra of Ib-c 
SNe, as was shown in Branch et al. 2002, Baron 
et. al. 2005, Branch et al. 2006, Elmhamdi et 
al. 2006, Parrent et al. 2007. 

A competitive hypothesis for description of 
the spectral feature with minimum near 6IOOA 
may be its interpretation as an absorption com- 
ponent of the PCyg profile of the line CII 6580A 
(Branch et al. 2006). But if the whole extended 
(>/~3OOi?0) stellar-wind envelope was seen at 
first in X-rays and subsequently as a UV flare 
in the XRF 060218/SN 2006aj afterglow, then it 
is hydrogen that should be connected at least 
to that part of the relic envelope which did not 
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evolve and refers to the stellar- wind stage of evo- 
lution of the core-collapse progenitor star, and 
which originated long before the SN explosion. 
The fact that it is hydrogen that always shows 
the largest contrast of velocity between the HI 
layer and the photosphere in comparison with 
other elements in spectra of Type Ib-c core- 
sollapse or "stripped-envelope" SNe (Branch et 
al. 2002 (see Fig. 23), Branch et al. 2001 (. Fig. 
9), Elmhamdi et al. 2006 (see Fig. 5)) also is in 
favour of the idea that it is these layers related 
to the wind envelope that are the first to come 
into motion result of SN explosion. One 
can estimate also the mass of this moving part 
of the envelope (the mass of gas in the Ha or 
HI layer) with the help of the equation from the 
paper Elmhamdi et al. (2006) in the following 
way: 

M{Mq) ~ (2.38 X W-^)v^h/T{Ha) , 
where time after the explosion t^^ is in days, 
the velocity of layer V4 is in units of 10^ kms~^, 
and T{Ha) is the depth of line formation. In 
this equation deduced for the Sobolev optical 
depth in the expanding envelope (Castor 1970; 
Elmhamdi et al. 2006), we used the velocity 
24000 km s^^ which refers only to that part of 
the HI layer which came into motion (see data 
in Table 3). This corresponds to the detached 
case in Fig. HI As a result, the mass of the HI 
layer turns out to be ~O.OOO6M0, and its dis- 
tance to the center by this time — 3.55 days 
after the beginning of the SN explosion — is in 
any case not less than 7.36 x lO^^cm. 



4.3 The Colgate shock-breakout effect 
in XRF060218/SN2006aj and in other 
SNe — hght curves, spectra, luminosities 
and sizes. 

The Type lb and Ic SNe have long been ob- 
served, it is has already long been considered 
that their progenitors are most probably the WR 
stars surrounded by more or less dense wind en- 
velope which resulted from evolution of a core- 
collapse star. The shock arising in explosion of 
evolved star core passes through the envelope 
and leads to a bright and short X-ray and UV 
flash which can last several hours — duration of 
the flare depends on how massive and extended 
was the wind envelope surrounding the progeni- 
tor star before the SN explosion. Interaction be- 
tween the shock caused by the SN explosion and 
the envelope (the shock breakout effect) has also 
been predicted and comprehended long ago (Col- 
gate 1968, Bisnovatyi-Kogan et al. 1975, Blin- 
nikov et al., 2002). In particular, the paper by 
Calzavara and Matzner (2004) is dedicated to fu- 
ture systematic observations of this effect. But 
so far, before the event XRF 060218/SN 2006aj, 
observations of the almost total effect has been 
gained for only small amount of the core-collapse 
SNe. 

The shock breakout effect which was very 
short due to the compactness of the blue su- 
pergiant (2O-3Oi?0) was observed (though not 
from the very beginning) in the famous Type 
II SN 1987A: there are quite a few data and the 
physics is described rather well — see the large 
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review by Imshennik and Nadezhin (1988). Dur- 
ing a long time (2-3 weeks), one could observe 
the shock breakout effect in a very extended (~ 
10^5 cm) envelope of Type Iln SN 1994W (see, 
for example, in Chugai et al. (2004)). In other 
cases of the Type Ib-c SNe the shock breakout 
effect was observed only in its very ending, be- 
fore the subsequent increasing of the SN bright- 
ness corresponding to the radioactive heating 
56Ni^56QQ^56pg _ fQj. SN1999ex, Stritzinger 

et al. (2002). Identical "remnants" of the 
shock breakout effect were seen (in the R band) 
even in SN 1998bw, which is usually related to 
GRB 980425 (see Galama et al, 1998). 

SN1993J in M81: It can be said that 
SN 1993J was observed almost in the very be- 
ginning of explosion. "Almost" because in this 
case the moment of the SN explosion beginning 
is known with an accuracy of ~12 hours, but not 
to seconds as in XRF/GRB SNe. Initially, this 
SN was classified by spectrum as the Type II SN 
because of appearance of hydrogen lines in its 
early spectra (Fig.[7|), but after a time, the SN 
type changed to lb (Flippenko, Matheson and 
Ho, 1993), when hydrogen in spectra stopped be- 
ing detected reliably. Early spectra of SN 1993J 
in Fig. [7| show a strong UV excess typical of the 
shock breakout effect and a smoothed comtin- 
uum almost without lines above 5000A, simi- 
lar to the spectrum of XRF 060218/SN 2006aj 
in Fig.m though in the case of SN2006aj the 
expansion velocities are much higher (see below 
the remarks about possible asymmetry of ex- 



plosion of SNe identified with GRBs). Besides, 
SN 1993J had an unusual light curve, when lu- 
minosity quickly increased to the first maximum 
and then headily fell during several days and 
then is has been slowly increasing for the sec- 
ond time during the subsequent two weeks. The 
behavior is analogous to that shown in Fig. [2] 
for XRF 060218/SN 2006aj. The SN1993J light 
curve has been long and well modeled by many 
groups: Nomoto et al. (1993), Young et al. 
(1995), Shigeyama et al. (1994). These pa- 
pers say that the shock breakout in SN 1993J is 
naturally explained by interaction between the 
shock and the extended (300/?© ~ 2.1 • 10^^ cm) 
hydrogen envelope of mass ~ IMq around the 
progenitor star. Here the luminosity achieved 
in the peak of the first/quick maximum which 
lasts only 4-5 hours can achieve a value of ~ 
-|^q45 gpg^g the total photon energy release 
~ 5 • 10^^ erg. 

No wonder that approximately iden- 
tical parameters of the envelope explain 
the shock breakout burst in the case of 
XRF/GRB 060218/SN 2006aj. Here also the 
total energy turns out to be of the same order 
of ~ 10^^ erg. But in the case of SN 1993J there 
was no gamma-ray burst, and, most probably, 
this is caused by asymmetry of the SN explosion 
(see below about that). 

If GRB/XRF 060218/SN 2006aj is another 
case when the Colgate shock breakout effect was 
observed in pure form from the very beginning 
of the SN explosion , then it may be understood 
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as quite a definite hint that the GRB itself is 
the first signal in gamma-rays that the collapse 
of a massive core started, which is followed by 
the total process — an anti-collapse and explo- 
sion of a massive SN: after the GRB an X-ray 
flash (XRF) and then a powerful UV flash can 
be observed. Then the key moment in opti- 
cal opbservations of transient sources related to 
XRF/GRBs can be the search for all manifes- 
tations of wind envelopes around core-collapse 
progenitor stars in early spectra and in photom- 
etry of the XRF/GRB afterglows, because in- 
teraction of the SN shock and these envelopes is 
the very first event in optical after the beginning 
of collapse of the massive stellar core. Here one 
may expect observing new effects related, for ex- 
ample, to the same asymmetry of explosion. 

4.4 Asymmetry of the Type lb and Ic 
SNe explosions. 

The burst XRF 060218/SN 2006aj was a clas- 
sical XRF event indeed (Heisc ct al. 2001, Cam- 
pana et al., 2006). The fact that in the case 
of usual and nearby SNe the explosion does not 
begin with a GRB is naturally explained by an 
aymmetric, axial-symmetric or bipolar (with for- 
mation of jets) explosion of the core-collapse 
SNc. Now one of the most popular conceptions 
(sec references in the paper by Sodcrbcrg et al., 
2005) proceeds from the idea that in the case of 
flashes of the XRF type an observer is out of the 
beam in which the most gamma-ray radiation is 
concentrated for one reason or another. 

The farther is an observer from the SN ex- 



plosion axis, the more of X-ray radiation and 
the less gamma-ray quanta are in the spec- 
trum of the flash — GRBs transform to X-ray 
Rich GRBs (hke GRB 030329) and become X- 
ray Flashes (Sokolov et al. 2006). When ob- 
serving at an angle close to 90° to the SN ex- 
plosion axis, no GRB is seen; one observes only 
an XRF (X-ray Flash) and then a powerful UV 
flash caused by interaction in the shock and the 
envelope surrounding the pre-SN as was in the 
case of SN 1993J. 

Thus, if an SN is observed close to the explo- 
sion equator (and this situation is the most prob- 
able) and if there is a sufficiently dense stellar- 
wind envelope around a massive collapsing star 
core, then only the shock breakout effect is to 
be observed in X-rays and in optical. In that 
case the contribution of GRB afterglow into a 
light curve of a "usual" SN can be unnoticeable. 
One way or another, but it must be much less 
than for classical GRBs observed close to the SN 
explosion axis (the least probable situation). In 
this connection, Filippenko et al. (2006) noted 
recently that a substantially asymmetric explo- 
sion can be a genetic feature of core-collapse SNe 
of aZZ types, though it is not clear yet if the mech- 
anism generating the GRB is also responsible for 
the star explosion. 

V. CONCLUSION 

The paper gives additional arguments in 
favor of the stellar-wind origin of the shock 



16 



breakout effect detected previously from 
Swift/XRT/UVOT data (Campana et al., 2006) 
on XRF/GRB 060218. In our optical spectra 
of the XRF/GRB 060218 afterglow we detected 
features interpreted as hydrogen lines which 
were also observed in early ESO Lick, ESO VLT 
and NOT spectra (Plan et al. 2006, Sollerman 
et al. 2006). Hydrogen was detected in spectra 
of a GRB afterglow for the first time, which is 
a direct sign of a relic wind envelope around a 
core-collapse progenitor star. 

Results of modelling two BTA spectra ob- 
tained in 2.55 and 3.55 days after explo- 
sion of SN2006aj related to the X-ray flash 
XRF/GRB 060218 are represented. The spec- 
tra were modeled in the Sobolev approximation 
with the help of the SYNOW code (Branch et 
al. 2001; Elmhamdi et al. 2006). In these early 
spectra of the Type Ic SN2006aj we detected 
spectral features interpreted as: 

(1) the Ha PCyg profile for velocities of 
~33000kms~^ — a wide and almost unnotice- 
able deformation of continuum in the range of 
pa 5600 — 6600A for rest wavelengths (z=0) at 
the first epoch (2.55 days), and 

(2) a part of the Ha PCyg profile in absorp- 
tion blucshiftcd by 24000 km s^^ — a wide and 
low-contrast spectral feature at ~ 6IOOA (rest 
wavelength) at the second epoch (3.55 days). 

Evolution of the same spectral features can be 
traces also by spectra of SN 2006aj obtained with 
other telescopes (Plan et al. 2006, Sollerman et 
al. 2006). Such Ha lines can directly confirm 



the existence of a stellar-wind envelope which 
was already observed during XRF/GRB 060218 
itself as a powerful black-body component first 
in X-rays and then in the optical spectrum — 
the so-called shock breakout effect (Campana et 
al., 2006). Thus, taking into account early ob- 
servations carried out with BTA and other tele- 
scopes (Table 1), it may be said that we observe 
evolution of optical spectra of the core-collapse 
SN 2006aj — a transition from the phase of the 
Colgate shock breakout effect to spectra of the 
phase of brightness increase corresponding to the 
radioactive heating. 

Our identification of these features with the 
Ha hydrogen line in early spectra of the Type 
Ic SN 2006aj was also confirmed by spectral ob- 
servations of usual (non-identified with GRBs) 
Type Ic and lb core-collapse SNe and interpre- 
tation of their spectra (in which hydrogen was 
also found) with the SYNOW code (Branch et 
al. 2002; Branch et al. 2006; Elmhamdi et al. 
2006). 

If interpretation of the thermal component 
in the spectrum of GRB/XRF 060218 as inter- 
action between the SN shock and the wind en- 
velope around the SN2006aj/XRF 060218 pro- 
genitor star will be confirmed by observations 
of afterglow of other bursts, then it will give 
a new impulse to development of the theory 
of GRBs themselves and of the core-collapse 
SNe. The intermediate redshift GRB/SNe 
are observed relatively rarely (Chapman et 
al. 2007), but they are the most informative 
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events (such as XRF /GRB 060218/SN 2006aj or 
GRB030329/SN2003dh) from the point of view 
of comprehension of relation between GRBs and 
SNe. And the key moment of the study of these 
transient sources may be the search for manifes- 
tations of wind envelopes around core-collapse 
progenitor stars of GRBs both in early spectra 
and in the photometry of GRB afterglows. 
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TABLE I: Early spectra of the supernova XRF 060218/SN 2006aj obtained before Feb. 23 2006 
UT with different telescopes. The spectrum time is time after XRF/GRB 060218 Swift trigger: 
Tsp is in days after 2006 Feb. 18.149 UT. Only spectra with a high signal/noise ratio are given. 
Early spectra obtained by Modjaz et al. (2006) with the 1.5-meter telescope FLWO in 3.97 days 
after the burst were not included because of their low signal/noise ratio. 

Telescope Tsp and 2006 UT References 

MDM (2.4ni) 1.95 days (Feb.20.097) Mirabal et al.2006 

BTA (6ni) 2.55 days (Feb.20.70) FatkhuUin et al. 2006 

ESO VLT (8m) 2.89 days (Feb.21.041) Pian et al. 2006 

BTA (6ni) 3.55 days (Feb.21.70) FatkhuUin et al. 2006 

NOT (2.56m) 3.78 days (Feb.21.93) SoUerman et al. 2006 

ESO Lick (3m) 4.01 days (Feb. 22.159) Pian et al. 2006 

ESO VLT (8m) 4.876 days (Feb. 23.026) Pian et al. 2006 
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FIG. 1: Observational spectra of the XRF 060218/SN 2006aj afterglow obtained with the 
BTA (FatkhuUin et al. 2006). UT times and observational times after the beginning of the 
SN explosion are indicated. Emission lines of the host galaxy (z=0.0331) are also shown. 
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FIG. 2: Optical (U, B, V, R, I) and infrared (J) light curves of XRF 060218/SN 2006aj 
(Jelinek et al. 2007). The first maximum corresponds to the UV flash/burst — the shock 
breakout effect (see the text). Spectra from Table 1 refer to the transitional region near 
minimum of the curves. Arrows point times of BTA spectra after the beginning of SN 
explosion. To(days) corresponds to 2006 Feb. 18.149 UT. 



FIG. 3: Field of the optical tranisent. BTA optical observations started on Feb. 20.647, 
2006 UT (Ri 60 hours after the GRB event) P.A. = 3° is a position angle of sht, 

VoT = 18.16, {B - R)oT = 0.3 




FIG. 4: Line profiles corresponding to the cases when envelope layers (i.e. layers over the 

photosphere) detach or do not detach from the expanding photosphere, when the gas 
expansion velocity increases proportionally to distance to the center {v ^ r, see the text). 
The shaded regions form the absorption component of the PCyg profile. 
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TABLE II: The table of model computation parameters in the case of layers undetached from 
the photosphere for all elements and their ions ("the undetached case" in Fig.|3|). The parameters 
correspond to the synthetic spectrum represented by the thick black line in Fig. [5] (see the text). 
More detailed information about setting parameters for the SYNOW code see in 
http: / / www.nhn.ou.edu/~parrent / synow.htmll 



Parameters 






"the undetached case" 


"Vphot " 






33000 


V ^ 
V max 






70000 








9000 


ai 


HI, 


He I, Fe HI, 


Fe II, Si II, Si I, I, Ca II, Ti II, CII, CI, NI, Mgll, Mgl, Mnll, Nel 


taul " 


•2, 


.08, .04, .04, 


.0002, .0002, .002, .5, .02, .0005, .0005, .1, .0002, . 0002, .0002, .0005 


vmine ^ 






33.00 for all ions 


ve 9 






20.00 for all ions 



''Vphot is velocity of the photosphere in km/s. 
^Vmax is the upper hmit of velocities in the model. 
'^tbb is temperature of the photosphere in Kelvin degrees, 
"^ai are ionization stages of all ions considered in the 

model. 

"^taul are corresponding optical depths of line formation 
of every ion. 

^vmine are the least velocities over the photosphere for 

every ion (in 1000 km s~^) 
*ve are the characteristic velocityes (ve) in the used law 

T(r) ~ exp(—v{r)/v^)) of relation between optical depth 

of lines of a given ion r and v, where v ^ r. 
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FIG. 5: The spectrum of the XRF/GRB060218/SN2006aj afterglow in rest wavelengths 
(z=0) obtained with BTA in 2.55 days and corrected for galactic extinction. Its fitting for 
"the undetached case" by synthetic spectra with the velocity of the photosphere (Vphot), all 

elements and their ions equal to 33000 km s~^ is shown by smooth lines differing only in 
the blue range of the spectrum at A < 4000^ (see the text). Main parameters of calculation 
of the synthetic spectrum represented by the thick black line are given in Table II. HI 
denotes the Ha PCyg profile at Vphot = 33000 km s^^. The model spectrum for the 
photosphere velocity 8000 km is shown for example by the dashed line as an example of 

the Ha PCyg profile. 
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TABLE III: Parameters of the SYNOW model for the case of layers of some elements 
undetached from the photosphere in the expanding SN envelope ("the detached case" in Fig.U]). 
The parameters correspond to one of the synthetic spectra shown in Fig. [6] by the thick black line 
(see the text). 



Parameters 


"the detached case" 




18000 


^ max 


75000 




8200 


ai 


HI, Hel, Fell, Sill, 01, Cal, Call, Till, CII, CI, NI, Mgll, Mgl, Nal 


taul 


0.16, 0.02, 0.08, 0.03, 0.1, 0.3, 0.3, 0.02, 0.0004, 0.04, 0.100, 0.005, 0.03, 0.02 


vmine 


24.00, 24.00, 24.00, 24.00, 18.00, 18.00, 18.0, 18.00, 24.0, 24.0, 18.0, 24.0, 24.0, 18.00 


ve 


10.00, 20.00, 20.00, 20.00, 20.00, 10.00, 10.0, 20.00, 10.0, 10.0, 20.0, 20.0, 20.0, 20.00 
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FIG. 6: The spectrum of SN2006aj/XRF060218 (rest wavelength at z=0) obtained with 
BTA in 3.55 days and corrected for galactic extinction. Synthetic spectra are shown by 
smooth lines. Locations of spectral lines of some ions and blends of their lines are shown in 
those parts of the spectrum where contribution of this ion into the spectrum is essential for 
given model parameters. The thick black line is the synthetic spectrum with parameters 

from Table III at which the absorption with minimum about 6100^ is described by 
suppressing influence of HI for "the detached case" . This is a strongly blue-shifted part of 
the Ha PCyg profile at the velocity of expansion of the detached HI layer equal to 

24000 kms-'^. 
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FIG. 7: The earliest optical spectrum of SN1993J (see the text). The spectrum is taken 
from the database SUSPECT - The Online Supernova Spectrum Database 
http://bruford.nhn.ou.edu/suspect/ (Richardson et. al. 2002) 



